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^ , Abstract 

^ ■ In our previous work of Ref. [HI we studied the stability of the RS2-model with a 

p\J , nonminimany coupled bulk scalar field (j), and we found that in appropriate regions of ^ 

^ I the standard RS2-vacuum becomes unstable. The question that arises is whether there 

■ exist other new static stable solutions where the system can relax. In this work, by 

. solving numerically the Einstein equations with the appropriate boundary conditions 

I on the brane, we find that depending on the value of the nonminimal coupling ^, this 

\^ • model possesses three classes of new static solutions with different characteristics. We 

. also examine what happens when the fine tuning of the RS2-model is violated, and we 

^ I obtain that these three classes of solutions are preserved in appropriate regions of the 

"V* ' parameter space of the problem. The stability properties and possible physical impli- 
cations of these new solutions are discussed in the main part of this paper. Especially 

r~| ! in the case where £, = £,c (Cc is the five dimensional conformal coupling) and the fine 

I tuning is violated, we obtain a physically interesting static stable solution. 
• ^ , 

X: 

c5 ; 1 Introduction 

A possible extension of conventional four-dimensional field theory models is based on the 
well known brane world scenario. According to this scenario, ordinary matter is assumed to 
be trapped in a submanifold with three spatial dimensions (brane world) that is embedded 
in a multi-dimensional manifold (bulk). Contrary to ordinary matter, gravitons are allowed 
to propagate in the bulk. The new feature in this scenario is that the extra dimensions can 
be large or even infinite. In addition brane world models predict interesting phenomenology 
even at TeV scale PP, and put on a new basis fundamental problems such as the hierarchy 
and the cosmological constant problem (see also the reviews and references therein). 

The typical example of a Brane-world model with an infinite warped extra dimension 
is the second Randall-Sundrum model (RS2-model) 0. In this scenario, we have a single 
brane with a positive energy density (the tension a), whereas the bulk has a negative five 
dimensional cosmological constant A. The corresponding Einstein equations have a solution 
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only if a fine-tuning condition is satisfied (A = — in units where SttG^ = 1. This solution 
implies a space-time geometry of the form of AdS^ around the brane, and Minkowski with 
zero effective cosmological constant on the brane. An extension of the RS2-model with a 
second negative tension —a brane, is the RSI-model j3]. In this case we have an orbifolded 
extra dimension of radius r^. The two branes are fitted to the fixed points of the orbifold, 
z = and Zc = vrr^ with tensions a and —a correspondingly. The particles of the standard 
model are assumed to be trapped on the negative tension brane, which is called visible, while 
the positive tension brane is called hidden. 

In this work we study the RS2-model with a nonminimally coupled bulk scalar field, via an 
interaction term of the form — |,^i?</)^, where ^ is a dimensionless coupling ^ The motivation 
is to examine whether this model possesses static solutions other than the standard one of 
the RS2- model (see Eq. (2) below). 

In particular we show, by solving numerically the Einstein equations with the appropriate 
boundary conditions on the brane, that according to the value of the nonminimal coupling ^ 
our model possesses three classes of new static solutions: (a) for ^ < the solutions develop 
a naked singularity in the bulk, (b) for ^ > (where C,c is the five dimensional conformal 
coupling) we find that the warp factor a{z) is of the order of unity near the brane and 
increases exponentially {a{z) ~ e^^), as 2; ^ +00, while the scalar field (j){z) is nonzero on 
the brane and tends rapidly to zero in the bulk (in this case the space-time is asymptotically 
AdS^), and (c) for < ,^ < the warp factor a{z) and the scalar field (l){z) tend rapidly to 
infinity. Contrary to case (b), where the scalar curvature is asymptotically constant, in case 
(c) the scalar curvature tends to infinity. In addition we examine what happens when the 
fine tuning of the RS2-model is violated, and we obtain that in appropriate regions of the 
parameters of the problem the three classes of solutions we described above are preserved. 

This work is motivated by our previous paper In particular in |2] we investigated the 
spectrum of a nonminimally coupled bulk scalar field in the background of the RS2-metric. 
We obtain that for ^ < the spectrum of the scalar field exhibits a unique bound state with 
negative energy, or a tachyon mode. Note that the existence of a localized tachyon mode 
indicates a gravity-induced Dvali-Shifman Mechanism jH] as we argue in It is worth to 
note, that this mechanism for the localization of Gauge field, from the point of view of lattice, 
has been investigated in Refs. [3 |H] . In |5j we had not examined the spectrum of the scalar 
field in the case of ^ > 0. We complete this investigation here. In the region < ^ < we 
obtain that there is no tachyon mode (see Appendix A). However the tachyon mode returns 
^ for ^ > ^c- The existence of tachyon modes for ^ < or ^ > implies an instability for 
both the RS2-metric and the scalar field (0 = 0) (see also 0). The question that arises 
is whether there exist other new static stable solutions where the system can relax. As we 
discussed in the previous paragraph this model indeed possesses static solutions different 
from the standard one of the RS2-model. The stability and possible implications of these 
numerical solutions are discussed in the main part of this paper. 

The spectrum of a nonminimally coupled bulk scalar field in the case of RSl-model j3] 
has been also considered [HlEl; where it was found that the tachyon character of the model 

^The coupling ^ possesses two characteristic values: a) the minimal coupling for ^ = and b) the 
conformal coupling for = 3/16 

^In this region of ^ the spectrum of the scalar field exhibits at least one tachyon mode, while for larger 
values of ^ it is possible to have two or more tachyon modes, see Appendix A 
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remains in appropriate regions of ^. However, the authors of Refs. [HI El are not interested 
to examine the stabihty of the RS-model, or to solve the Einstein equations. They mainly 
use this result in order to put the standard model on the brane (negative tension brane) 
with a bulk Higgs field and a gravity-induced Higgs mechanism. 

In the case of four dimensions the same nonminimally coupled model has been considered 
in Refs. [TTl [T2] in the background of de Sitter space-time. It is obtained that for a specific 
range of values of ^ the scalar field is rendered unstable, and this result has straightforward 
implications to the cosmological constant problem. 

2 RS2-model with a nonminimally coupled scalar field 

The action which describes the RS2-model, if we set SttG^ = 1 {G^ is the five-dimensional 
Newton constant), is 

S = Jd'xCns = J d'x Q^(i?-2A)-a5(^)v/|^7(^— )|) (1) 

where d^x = d'^xdz, and z parameterizes the extra dimension. In addition R is the five- 
dimensional Ricci scalar, g is the determinant of the five- dimensional metric tensor qmn 
{M,N = 0,1,. ..,4), and 

gbrane jg determinant of the induced metric on the brane. We 
adopt the mostly plus sign convention for the metric ISj. 

If the fine tuning A = is satisfied, the Einstein equations have a stable static solution 
of the form 

ds'^ = e~'^^^'''\-dxl + dx\ + dxl + dxl) + dz^ (2) 



where k = 




In this work we aim at the studying of the RS2-model with a nonminimally coupled bulk 
scalar field. The action of this model is: 

S = j d^'x {Crsf + C^) (3) 

The gravity part of the lagrangian is given by the equation 

Crsf = {F{ct>)R - A) - a5{z)^W^\ (4) 

where the factor 

m = \ii-i^') (5) 

corresponds to a nonminimally coupled scalar field with an interaction term of the form 
i^int = — l^-R 0^, and ^ is a dimensionless coupling. 
The scalar field part of the lagrangian is 

= ^(-^/'^Vm0V;v0-^(0)) (6) 
where the potential is assumed to be of the standard form V"(0) = Xcj)"^. 
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The Einstein equations, which correspond to the action of Eq. (3) are 



^g{brane) | 

GMN + ^gMN + cr6{z)^ — 1= — QixAS"^ = T^Slf (7) 

V l^'l 

where the energy momentum tensor for the scalar field is 

T^M^ = VM0ViV0-^MAr[^/''Vp0Vs0+y(0)]+2VMVA.F(0)-25MAfnF(0) + (l-2F(0))G'M7V 

(8) 

The equation of motion for the scalar field is 

, dFU)^ dVU) ^ 

The above equation is not independent of the Einstein equations (7), as it is equivalent to 
the conservation equation V^Tj^lf = 0, where T^]^ is given by Eq. (8). 
We are looking for static solutions of the form 

ds'^ ^a'^(z)(-dxl + dxl + dxl + dxl)+dz'^, (f) ^ (f)(z) (10) 

From Einstein Equations (Eq. (7)) we obtain two independent equations: 

Gii + Agu + a6{z)gu^T^f^ (11) 
G,, + Ag,,^T^P (12) 

where i = 0, 1, 2, 3. If we set a{z) = e^^^^ and use Eqs. (8), (10), (11) and (12) we get 

Fi = : 3(1 - ^d^'iz)) {A"iz) + 2A'{zf) + A + (^ - 2O0'(^)' + ^(0(^)) 

-2C<P{z)<P"{z) - 6iA'{z)<P{z)<P'{z) + a6{z) = (13) 



F2 = : 6(1 - C<f>\z))A'{zr + A - ^<P'{z)' + V{<t>{z)) - 8CA'iz)<Piz)<f>'iz) = (14) 

From Eq. (9) for the scalar field we get 

F3 = : -(l>"(z) - AA'(z)(l)'(z) - e {8A"(z) + 20A'(z)^) (l>(z) + V'(^) = (15) 

We can show that Eq. (15) (or equation F3 = 0) can be found from Eqs. (13) and (14) (or 
equations Fi = and F2 = 0) by taking the combination —4A'{z){Fi — F2) + -Fg = 0. Note 
that the relation -AA'{z){Fi - F2) + = is equivalent with the condition V^Tj^lf = 0. 

As we have already mentioned, Eqs. (13), (14) and (15) are not independent. The solu- 
tions can be obtained by integrating the second order differential equations (13) and (15). 
The first order differential equation (14) is an integral of motion of Eq. (13) and (15), and 
acts as a constraint between A{z), (j){z) and their first derivatives. 
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In particular we choose to solve the second order differential equation Fi — F2 = (see 
Eqs. (13) and (14)) and F3 = (see Eq. (15)), or 



3(1 - ^cf>'{z))A"{z) + (1 - 20cl>'{zf - 2^cf>{z)cf>"{z) + 2^A'{z)cP{z)cf>'{z) + aS{z) = (16) 

- (j)"{z) - 4A'{z)(l)'{z) - e {8A"{z) + 20A'{zy) (j){z) + V'{^) = (17) 

As this system is complicated we will not look for analytical solutions, but we will try to 
solve it numerically. For the numerical integration of Eqs. (16) and (17) it is necessary 
to know the values of ^(0), 0(0), ^'(O''") and 0'(O''"). These values are determined by the 
junction conditions (see Eqs. (18) and (19) below) and the constraint of Eq. (14). 

The delta function in Eq. (16) implies that the first derivatives of A{z) and (f){z) are dis- 
continuous on the brane (z=0). If we integrate Eqs. (16) and (17) over z, in an infinitesimal 
interval [— e, e], we get the junction conditions: 

3(1 - e0'(O)) {A'{0+) - A'{0-)) - 2^m (0'(O+) - 0'(O-)) + <7 = (18) 

- (0'(O+) - 0'(O-)) - 8e (^'(0+) - A'(0-)) m = (19) 

If we take into account the Z2 symmetry we found A'{0^) — A'{0~) = 2A'[0'^) and 0'(O"'") — 
(j)'{0') = 20'(O+), thus 

6(1 - e0'(O))^'(O+) - 4e0(O)0'(O+) + (7 = (20) 
0'(O+) + S^A'{0^)m = (21) 

By solving these equations with respect to A'(0+) and to 0'(O^) we obtain 

^ " (6 - 6^0(0)2 + 32^0(0)2) ^^^^ 

6'm = 8ea0(O) 

' (6-6^0(0)2 + 32^20(0)2) ^ ' 

However, the problem has an additional constraint for the first derivatives of A(z) and (f){z), 
which is given by the first order differential equation (14) (or F2 = 0). Thus if we replace 
Eqs. (22), (23) in Eq. (14) we obtain the following sixth order algebraic equation for 0(0): 



a2 



(6-6^0(0)2 + 32^20(0)2) 



A + V(0(O)) = O (24) 



Note that if we set 0(0)2 ~ x in Eq. (24) we obtain a third order algebraic equation (if 
1/(0) = A0^). If additionally the fine tuning A = is satisfied, the constant term of 
the third order algebraic equation (24) is zero, and thus we have a second order algebraic 
equation to solve, in order to determine the value of 0(^) on the brane. 
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Figure 1: The warp factor a{z) and the scalar field (f){z), in case (a) < 0), as a function 
of kz for a/k = 6, A/fc^ = —6, ^ = —0.1, X/k"^ = 0.01. The values for a, A satisfy the fine 
tuning of the RS-model. We see that the metric of our model exhibits a naked singularity, 
as the warp factor vanishes for kzs — 1.197. 

3 Numerical results 

For the determination of the two unknown functions A{z) and (j){z) we can solve the system 
of second order differential equations (16) and (17) for z > numerically. In order to 
integrate it is necessary to know the values of A{0), 0(0), ^'(0"*") and 0'(O+). If we assume 
that the warp factor is normalized to unity on the brane (or a(0) = 1) wc find that A{0) = 
(note that a{z) = e^'^^^). The value of 0(0) is obtained by solving Eq. (24), and the values 
of A'{0^) and 0'(O+) can be found from Eqs. (22) and (23). Then it is an easy task to 
use a routine of Fortran or Mathematica to extract the numerical solutions. Note that the 
numerical results we obtain satisfy with a very good accuracy also the first order differential 
equation (14). 

The model we examine has four independent parameters ^,X,A,a. We will keep fixed 
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the parameters A, A, a, assuming that the fine tuning A = — ^ is satisfied, and we will vary 
the parameter ^. As we discuss in the following sections, depending on the value of ^ we 
find three classes of numerical solutions with different characteristics. Also we investigate 
what happens when the fine tuning is violated, and wc find that in appropriate regions of 
the parameter space, the three classes of solutions wc described are preserved. However, 
there are regions of the parameters where there are no static solutions (or Eq. (24) has no 
real solutions) . A thorough investigation of these regions is very extended and it is beyond 
the scope of this paper. 



It is convenient in numerical analysis to perform the rescaling y — kz where k — 
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Figure 2: The warp factor a{z) and the scalar field (f){z), in case (a) < 0), as a function 
of kz for a/k = 6, A/k^ = -6, ^ = -0.001, X/k^ = 0.01. The values for a, A satisfy the fine 
tuning of the RS-model. We see that the metric of our model exhibits a naked singularity, 
as the warp factor vanishes for kzg = 1.994. 

Then the Einstein equations (13), (14) remain unchanged if the parameters A and a are 
divided by fc^ and k correspondingly. Note that A/k"^ = —6. The equation (15) for the 
scalar field remains unchanged if the parameter A is divided by fc^. 

We emphasize that the system of Eqs. (16) and (17) with the boundary conditions of Eq. 

2 

(22), (23) and (24), in the case of the fine-tuning A = — has an obvious analytic solution 
of the form of Eq. (10) with an exponential warp factor a{z) = e"'^'^' and scalar field vacuum 
equal to zero (j){z) = 0. This solution is identical to the well known solution of Eq. (2) for 
the RS2- model. However, it is unstable for ^ < and C, > (see Ref. and Appendix A), 
thus it is worth investigating whether this model posses static solutions other than that of 
Eq. (2), as we do in the rest of this section. 

3.1 Case (a) < 0) 

The main feature of the numerical solutions for ,^ < is a naked singularity at finite proper 
distance Zs in the bulk. The scalar field 0(z) is almost constant near the brane, and tends 
to infinity as z tends to the singularity point in the bulk. 

In Fig. ^we have plotted the warp factor a{z) and the scalar field (plz) as a function of z 
for a/k = 6, A/fc^ = —6, ^ = —0.1, X/k'^ = 0.01. The values for a, A satisfy the fine tuning 
of the RS-model. We see that the metric of our model exhibits a naked singularity, as the 
warp factor vanishes for kzg = 1.197. Note that the warp factor for ^ = —0.1 is completely 
different from the exponential profile of the warp factor of the RS2-model. 
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Figure 3: The warp factor a{z), in case (a) (.^ < 0), as a function of kz for two values of 
a/k = 6.5 and 5.5, A/k"^ = —6, C, = —0.1, X/k"^ = 0.01. The values for a, A do not satisfy 
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the fine tuning of the RS-model. Note that even when the fine tuning A = — ^ is violated 
the induced metric on the brane is Minkowski. 

However for very small absolute values of ^ (for example ^ = —0.001), as we see in Fig. 
the warp factor near the brane is almost identical with the exponential profile of Eq. (2) 
{a{z) = e~^^), (in Fig. |21we we have plotted the warp factor a{z) and the scalar field (f){z) 
as a function of z for a/k = 6, A/k"^ = —6, ^ = —0.001, X/k"^ = 0.01 and k = 1). However 
even in this case we can not avoid a naked singularity in the bulk at Zsk = 1.994. We have 
checked that a naked singularity in the bulk appears for arbitrarily small negative values of 

In Fig. El we examine the case where we have a violation of the fine tuning A 7^ — As 
a consequence, the RS2-metric of Eq. (2) (exponential warp factor) is not a solution of the 
RS2-model with a nonminimally coupled bulk scalar field. However, we see (Fig. ^ that 
this model possesses static solutions of the form of Eq. (10) with a naked singularity in the 
bulk. It is interesting to note that the induced metric on the brane is Minkowski, even when 

2 

the fine tuning A = — ^ is violated. We remind the reader that in the case of a violation 
of the fine tuning, the RS2-model (without the nonminimally coupled scalar field) possesses 
solutions with an AdS^ or dS^ induced metric on the brane, see for example Ref. Of 
course if we wish to construct a brane world scenario, we should put a second brane in the 
bulk before the singularity. However, in this case the tension of the second brane a', and 
the parameters A and a of the first brane should be finely tuned suitably, if we want the 
Einstein equations and the boundary conditions on the second brane to be satisfied (see for 
example Ref. j2] and references therein). 
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Note that the numerical solutions we found for ^ < are similar to the analytic solutions 
in Ref. [T3]. However, the model in Ref. is quite different from that we assume here. 
In Ref. [13] it is argued that this kind of solutions may have physical interest even without 
a second brane. In particular, there is a possibility that quantum gravity effects near the 
singularity push the singularity to infinite proper distance in the bulk. In this way, the 
singularity acts effectively as a physical end to the extra dimension z (for more details see 
Ref. IIS]). A way to take into account quantum gravity effects, is to add to the gravity 
action of Eq. (1) a Gauss-Bonnet term. We have performed numerical computations also for 
this case and the answer is that Gauss-Bonnet gravity can not resolve the naked singularity 
(see also jTHI)- The naked singularity remains for arbitrary values (positive or negative) of 
the free parameter in front of the Gauss-Bonnet term. 

In our previous work 5j we pointed out the instability of both the RS2-metric of Eq. 
(2) and the scalar field vacuum = 0, for ^ < 0. Responsible for this instability is the 
existence of a unique tachyon mode localized on the brane. By using this result in Ref j3] 
we tried to guess the form of the new static stable solution where the system is expected to 
relax. In particular, we supposed that the profile of the scalar field vacuum in the bulk is 
proportional to the wave function of the tachyon mode (or the scalar field vacuum is nonzero 
on the brane and tends rapidly to zero in the bulk). If we take for granted that (f){z) —>■ for 
1^1 +00, the warp factor for \z\ +oo must be of the form a{z) = e"*^'^' or a{z) = e*^'^', 
as only these two functions satisfy the Einstein equations if we set 0(z) = 0. In Ref. [HI we 
guess the most plausible behavior for the new static stable metric. However in this work, 
by solving numerically the Einstein equations with the appropriate boundary conditions, we 
obtain a completely different behavior for the warp factor and the scalar field vacuum. The 
scalar field vacuum has not the profile of the tachyon mode. As we see in Fig. ^ and Fig. |21 
the warp factor vanishes in finite proper distance in the bulk, creating a naked singularity, 
whereas the scalar field vacuum tends to infinity as we approach the naked singularity. 

Another way to try to resolve the naked singularity is to include a suitable energy mo- 
mentum tensor Tmn in the right hand side of Einstein equations (7). We have observed 
that an energy momentum tensor of the form Tzz = c/a{z)'^, with all the other components 
equal to zero, is possible to give stable solutions of the form we described in the previous 
paragraph. However, we could not find a physical explanation for an energy momentum 
tensor of this form. 

3.2 Case (b) > Q 

For ^ > "we obtain a different class of numerical solutions. As we see in Fig. El the warp 
factor a{z) is of the order of unity in a small region near the brane and increases exponentially 
{a{z) —>■ e^^), as z +oo. From this figure we can estimate that the warp factor a{z) is 
of the order of unity for < 2; < 3.5 and the pure exponential behavior a{z) = e^^ begins 
for z > b. The exponential behavior of the warp factor is confirmed with great accuracy by 
numerical computations, as in the left-hand panel of Fig. El we observe that A'{z) ~ 1 for 
large values of z. The scalar field (l){z), as we see in Fig. El is nonzero on the brane (in a 
rather small region (0 < < 4), and for kz > 5 tends rapidly to zero. We have checked 
carefully that this class of solutions appears even when ^ exceeds by a small amount the 
conformal coupling (^ = + SC. with SC, « 1). 
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Figure 4: The warp factor a{z), in case (b) > C,c), as a function of z for two values of 
a/k = 6 and 2, A/A;^ = —6, ^ = ^c + ^, ^/k'^ = 0.01. Note that the solution exists even when 
the values of a and A do not satisfy the fine tuning of the RS-model. 

2 

Note that this class of numerical solutions is preserved even when the fine tuning A = — ^ 
is not satisfied (see the down panel of Fig. [Hand the right panel of Fig. = 2, A = —6)). 
In this case also, A'{z) ~ 1 for large values of z. We have checked numerically that for fixed 
negative ^ A there are solutions of the form of Figs. ID and IHl only if the brane tension a is 
smaller than the absolute value of the five dimensional cosmological constant A. 

A feature of this class of solutions is that the warp factor tends to infinity a{z) — ^ e'^^ for 
large z. If we wish to construct a brane world scenario, it is necessary to include a second 
positive tension brane far away from the first, where the numerical value of the scalar field 
is practically zero. In this case the tension a' of the second brane should satisfy the fine 
tuning condition A = — |— . 

However the most important topic for the construction of a realistic brane world scenario 
is the stability of this class of solutions. In order to study the stability, we will replace 
— > + in the lagrangian 

^ = ^(-^^?^^Vm0V;v0-^W-V^(0)) (25) 

where is a small perturbation around the classic solution. If we keep only second order 
terms of 0, we obtain that an effective mass ^e// ~ ^-^ + 12A0(2;)^ arises for the scalar 
field perturbation 0. In the right-hand panel of Fig. El we have plotted the effective mass 
•^Note that for nonnegative values of A there are no solutions of the form we described. 
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Figure 5: The A'{z)/k and mlj-jr/k'^, in case (b) > ^c), as a function of kz for a/k = 6, 
A//c^ = —6, ^ = + 1, A//c^ = 0.01. Note that asymptotically both the derivative A'{z) and 
the miff tend to a constant value, or the space-time is asymptotically ^^5*5. 



^eff ~ + 12A0(2;)^ as a function of z. We see that as 2; ^ +00 the square of the effective 
mass tends to a constant negative value. This result indicates ^ an instability of the solutions 
of class (b) against scalar field perturbation. Hence this class of solutions can not be used 
for the construction of realistic brane world scenarios. 



3.3 Case (c) (0 < <^ < ^c) 

For < ^ < the standard RS2-metric is stable, as there are no tachyon modes for the 
scalar field perturbations around this solution (see Appendix A). As we have mentioned, in 
this case a second static solution can be found by solving the Einstein equations numerically. 
This second solution exists only if Eq. (24) has at least a nonzero real root for 0(0). As we 
see in Fig. [7|for < ^ < the warp factor a{z) of this solution tends rapidly to infinity 
(faster than case (b)), also the scalar field (j){z) is nonzero on the brane and tends rapidly 
to infinity in the bulk. Contrary to case (b), where the space-time is asymptotically AdS^, 
in this case the scalar curvature tends to infinity. Solutions of this kind exist even when the 
fine tuning of the RS-model is violated. 

If we wish to construct a brane world model, we should put a second brane in the bulk. 
However, in order to satisfy simultaneously the three boundary conditions of Eqs. (22), (23) 

''A negative mass term is not enough to establish the instabihty of the model. In the case of curved 
space-time it is necessary for at least one tachyon mode for the scalar field to exist. In Appendix B we show 
that indeed the spectrum of the scalar field exhibits a tachyon character. In particular, there is a continuous 
spectrum of tachyon modes for ^ > ^c- 
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Figure 6: The scalar field (f){z), in case (b) > ^c), as a function of kz for two values of 
a/k — 6 and 2, A/A;^ = —6, ^ = + 1, — 0.01. Note that the solution exists even when 
the values of a and A do not satisfy the fine tuning of the RS-model. 

and (24) on the second brane, it is necessary to assume a fine tuning for the parameters 
of the model. Moreover, as we did in case (b), we have computed the effective mass of the 
scalar field around this solution and we have seen that it is positive. This implies that this 
class of solution is stable against perturbations of the scalar field. 

4 Minimal coupling and conformal coupling 

In this section we examine separately the special cases of ,^ = and ^ = ^c- Especially in the 
case of conformal coupling, when the fine tuning is violated, we have a physically interesting 
static stable solution. 

4.1 Minimal coupling ^ = 

The algebraic equation (24), for 0(0), can be solved analytically for ,^ = or ,^ = If we 
set ^ = or ,^ = we find the equation A0(O)*^ + a^/6 + A = 0. In the case of fine tuning we 
obtain 0(0) = 0, hence the solution of the model we examine is the standard one with warp 
factor a{z) — e~^\^\ and ^(z) = 0. If (7^/6 + A > 0, the above equation has no real roots and 
our model has no static solutions. On the other hand, for cr^/ 6 + A < we have two real 
roots 0(0) = ±-y/— A~^((j2/6 + A), and in this case the static numerical solutions for ^ = 
have the same characteristics with those of class (a) (^ < 0) with a naked singularity. 



12 




Figure 7: The warp factor a{z) and the scalar field (j){z), in case (c) (0 < ^ < ^c), as a 
function of kz for a = 6, A = —6, C, = 0.1, A = 0.01. The values for a, A satisfy the fine 
tuning of the RS-model. 

4.2 Conformal coupling ^ = 

2 

In this Section we examine the model when ^ = and A 7^ — As we have already 
discussed in Section 4.1, if a/k < 6, the algebraic equation (24) for 0(0) has two real 

roots 0(0) = ±00, where 0o = ^ —X^^^a"^ /Q + A). In this case, the numerical solutions of 
the model exhibit the same characteristics with those of class (b) (^ > ^c)- In particular, 
the warp factor a{z) is of the order of unity near the brane and increases exponentially 
(a(z) ~ e'^'^' ), as z — > ±cxd, while the scalar field 0(z) is nonzero on the brane and tends 
rapidly to zero in the bulk, as we can see in the upper panels of Fig. |H| The exponential 
behavior of the warp factor a{z) = e^^^^ is confirmed in the down-left panel of Fig. |H1 The 
solutions in the case of conformal coupling are of the form we described only for suitable 
small values of A//c^ depending on the exact value of a (i.e. for a = 5 we have A/fc^ < 0.1). 
For larger values of A we find a different class of solutions, however it is not of physical 
interest . 

In the down-right panel of Fig. |H1 we observe that the effective mass of the scalar field 
is negative. However the solutions for ^ = are stable, as in Appendix C we find that 
the scalar field does not exhibit four dimensional tachyon modes. The result of stability 
implies that the solutions for = can be used as the background for the construction of 
a brane- world scenario. 

Now we will try to discuss possible physical implications of the solutions in the case of 
conformal coupling. As we have already mentioned, the warp factor behaves as e'^'^' for large 
z (see the down-left panel of Fig. IH)). It is well known that an exponentially increasing 
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Figure 8: a{z), 4>{z), A'{z)/k and mljj-{z)/k'^, for ^ = ^c, as a function of kz for a/k = 5, 
A/fc^ = —6 and A/A;^ = 0.01. Note that asymptotically both the derivative A'{z) and the 
mgjj tend to a constant value, or the space-time is asymptotically AdS^. 

warp factor comes from negative tension brane with —a'/k = —6. Hence, the positive 
tension brane ai z = {a/k = 5) together with the negative energy density of the scalar 
field vacuum Tgo^, act effectively as a negative tension brane, which creates the exponential 
profile e'^'^' for the warp factor. Note that the energy density of the scalar field has a finite 
size extension into the bulk of the order of 1/k, as we see in the upper-right panel of Fig. 
IHl The model is completed if we include a second positive tension brane with a'/k = 6, in 
a position Zc where the scalar field vacuum is practically zero [z^ » 5/k). This two-brane 
set up is very similar to the well known RSl-model (see introduction). In analogy with the 
RSl-model, we will assume that the standard model particles live in the effectively negative 
tension brane, or visible brane . Note that in our case the visible brane is a formation of the 
positive tension brane with a/k = 5 and the negative energy density Tqq ^ of the scalar field 
vacuum. 

The advantage of this model is that it incorporates a mechanism for the localization of 
standard model particles on the visible brane. We can enrich our model with a Gauge field 
symmetry, i.e. SU(5) (see also We assume that this gauge field symmetry (SU(5)) 

is spontaneously broken to SUc{3) x SUl{/2) x f/y(l) near the brane , due to the nonzero 
value of the scalar field, while it is restored in the bulk for z » where the value 
of the scalar field is practically zero. This phase structure, Higgs phase on the brane and 
confinement phase in the bulk, triggers the Dvali-Shifman mechanism j6| for localization of 



14 




Figure 9: a{z), 4>{z), A'{z)/k and TQQ{z)/k'^, for ^ = ^c, as a function of kz for a = 0, 
A/k^ = -6 and X/k"^ = 0.01. 

Gauge fields. In this way the gauge fields of SUc{S) x SUl{2) x f/y(l) and the matter fields 
with gauge charge (see Ref. j2]) are locahzed on the brane, as for escaping in the bulk, it 
requires energy equal to Agap, where Agap is the mass gap emerging from the nonperturbative 
confining dynamics of the SU(5) gauge field theory in the bulk. 

In the case of the RSl-model, when ordinary matter is localized on the negative tension 
brane, we have the relation Mp = ^{e^^""" - 1) (see Ref. 0), where Mp = lO^^Gel^ is 
the four-dimensional Planck scale, is the fundamental five- dimensional gravity scale, k 
is the inverse AdS^ radius, and Zc is the position of the second positive tension brane in the 
bulk. We have checked that the above relation is valid approximately also for the model we 
examine. If we wish to have strong gravity at TeV we must choose M^, = ITeV. According to 
this localization mechanism the particles on the visible brane see an effective length toward 
the extra dimension of the order of 5/k. If we take into account that the masses of ordinary 
particles are smaller than ITeV, we must choose k/5 » ITeV. For example, ii k = lOOTeV^ 
we obtain lOMp = M^e^^'', and thus the position of the second brane is Zc ~ 39/ k. 

Note that even for a zero tension {a = 0) brane at 2; = 0, we obtained solutions of the form 

we described in the previous paragraphs, see Fig. In this case the visible brane is formed 

only from the negative energy density of the nonminimally coupled bulk scalar field, as we 

see in the down-left panel of Fig. IHl In addition, we have found that the numerical solutions, 

for small values of X/k"^ (-^/^^ ^ 0.01), can be approximated with a very good accuracy 

2 i_ 

from the analytical expressions a(z) = fcosh(^A;|z|) j ^ and 0(z) = 0(0)/ fcosh(^A;|z|) j . 
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These expressions satisfy the constraint equation (14). However, their second derivatives are 
not in agreement with the second derivatives of the numerical solutions, hence the equations 
(16) and (17) are not satisfied by the above mentioned analytical expressions. 

5 Conclusions and discussion 

In this work we studied for the first time the RS2-model with a nonminimally coupled bulk 
scalar field, via an interaction term Cint = —\^R<p'^- By solving numerically the Einstein 
equations with the appropriate boundary conditions on the brane, we showed that depending 
on the value of the nonminimal coupling ^ this model possesses three classes of new static 
solutions with different characteristics. 

Class (a) < 0) develops a naked singularity in the bulk. Class (b) (^ > ^c) consists 
of solutions which are characterized by an exponential warp factor a{z) ~ e'^^, as z — >■ +cxd, 
and a scalar field (j){z) which is nonzero on the brane and tends rapidly to zero in the bulk. 
The solutions of class (c) (0 < ,^ < ^c) are characterized by a very fast increase of the warp 
factor (faster than class (b)), and a scalar field (j){z) which tends rapidly to infinity for large 
z. The cases of minimal and conformal coupling (^ = and ^ = ^c) have been discussed 
separately in Sections 4.1 and 4.2. 

An interesting point is that these three classes of solutions exist even when the standard 
fine tuning of the RS-model is violated. In order to construct a brane world model by using 
these solutions, it is necessary to assume a second brane in the bulk. However, in this case 
we should impose a new fine tuning between the parameters of the model, if we wish to 
satisfy the boundary conditions on the second brane. 

In addition, we have examined the stability properties of the new solutions. The solutions 
of class (a) {S, < 0) and class (b) {S, > C,c) are unstable against scalar field perturbations, as 
in Appendix B we have found that the spectrum of the scalar field around these solutions 
exhibits a tachyon character. The only way to render these solutions stable is to put a second 
brane in the bulk before the potential for the scalar field V(w) becomes negative (for details 
see Appendix B). The solutions of class (c) are stable, as the scalar field spectrum has no 
tachyon modes. 

We emphasize that this work has been motivated by our previous work of Ref . [3] , where 
we had assumed the same model with an additional gauge field symmetry for the scalar field. 
In particular, in that work we have argued for a gravity-induced localization mechanism for 
^ < 0, which is very useful for the localization of gauge fields on the brane. In this mechanism 
we have a specific phase structure: confinement phase in the bulk and Higgs phase on the 
brane. Gauge fields, and more generally fermions and bosons with gauge charge (see Refs. 
012), can not escape into the bulk unless we give them energy greater than the mass gap 
Agap, which emerges from the nonperturbative confining dynamics of the gauge field model in 
the bulk. In Ref. we observed that the effective mass mljr j = ^R, foi ^ < 0, is negative on 
the brane and positive in the bulk. This result indicates a phase structure which can trigger 
the Dvali-Shifman mechanism in a gravitational way. We would like to emphasize that the 
sign of the effective mass, in the case of curved space-time, is just an indication for the 
expected phase structure and is not a strict proof. A realistic gravity- induced localization 
mechanism requires a static stable solution with a nonzero scalar field on the brane which 
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vanishes rapidly in the bulk. In Ref. we assumed the existence of a solution of the form 
we described above. However, the correct way to find if a solution of this kind really exists, 
is to solve the complete system of Einstein equations, as we do in this paper. We found that 
no static stable solution of the required form exists for ^ < 0. Additionally, we obtained that 
for ^ < the solutions suffer from a naked singularity in the bulk, see Figs. I1I2I and El We 
tried to resolve this naked singularity by adding a Gauss-Bonnet term to the gravity action. 
We performed numerical computations and we found that the naked singularity remains 
even in the case of Gauss-Bonnet gravity. According to the above analysis, it seems that no 
realistic gravity- induced localization mechanism for ^ < can be constructed. 

Finally, we have examined the case of conformal coupling = ^c) when the fine tuning 
of the RS-model is violated, see Section 4.1. We obtain that the solutions for ^ = have the 
same characteristics with those of the solutions of class (b) < ^c)- However, there is an 
important difference between them, as the solutions for ^ = are stable against scalar field 
perturbation (see Appendix C), contrary to the case of second class, where the solutions are 
unstable (see Appendix B ). In Section 4.2 we argue that this class of static stable solutions 
can be used for the construction of a realistic brane world scenario, which is very similar to 
the well known RSl-model. In this scenario, the visible brane is a formation of the positive 
tension brane o" at 2; = 0, plus the negative energy density of the nonminimally coupled 
scalar field. The advantage of this model is that it incorporates a mechanism (gravity- 
induced Dvali-Shifman mechanism) for the localization of standard model particles on the 
visible brane together with the spontaneous breaking of a Grand Unified Gauge group. We 
would like to emphasize that the above mentioned scenario remains even in the case of a 
zero tension brane at 2; = (cr = 0). In this case, the visible brane is formatted only from 
the negative energy density of the bulk scalar field. 
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A Appendix: Spectrum of scalar field perturbations 
around the RS2-vacuum 

In Ref. jS| we studied the spectrum of scalar field perturbations in the background of the 
RS2- vacuum for C, < 0. In this appendix we complete this investigation by including positive 
values for ^. 

If we use the transformation 

((,k\z\ _ I) 

w = sgn{z)- (26) 
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the RS2-metric of Eq. (2) can be put into the manifestly conformal to the five- dimensional 
Minkowski space form 



where 



ds^ = a{wY{—dxl + dx\ + dx^ + dx\ + dw'^) 

1 



a{w) 



k\w\ + 1 



(27) 



(28) 



The lagrangian of the scalar field if we include an interaction term —\R(j)^ is written as 



L = ^J\9\[ -^5*'^Vm0V^0 - \iR <P' - V{<P) 



(29) 



If we consider a small perturbation around the scalar field vacuum (0 = 0) we find the 
corresponding linearized equation: 



1 



'M 



g\g^''d^<j>{x,w) 



We can set 



0(x, w) — e 



ipx 



+ CR{w) ${x,w) = 



(30) 



(31) 



where a{w) = l/(A;|-u;| + 1), and = Pi^P^ is the effective four dimensional mass. 
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Figure 11: For ^ = + 1 the spectrum of the scalar field contains exactly two tachyon modes 
with negative energies m\/k'^ = —0.278 and m\/k'^ = —1.163. In the figure we have plotted 
the wavefunctions ipi and ■?/'2 of the two tachyon modes as a function of z. 



The function tplw) satisfies the Schrondiger like equation 

- -^"{w) + \v{w) - m^] ip{w) = 



where the potential V{w) is equal to 

V{w) = ■ 
From Eqs. (27) and (32) we get 



a 



3/2, 



a 



3/2 1 



w 



+ ^a^{w)R{w) 



(32) 



(33) 



Viw) = -16ki^-^,){-Siw) + 



5k 



4:{k\w\ 



(34) 



where = 3/16 is the five dimensional conformal coupling. 

Note that the coefficient in front of the potential change sign when crosses the five 
dimensional conformal coupling. This result implies that the potential has two characteristic 
forms, as we see in the left-hand panel (^ < ^c) and the right-hand panel (^ > ^c) of Fig- CHI 

In the first case, where ^ < the coefficient of the delta function is negative. In Ref. 
[S] we have shown that if ^ < 0, the spectrum of the scalar field contains a unique tachyon 
mode localized on the brane (see Eq. (32) in Ref. ^). It is well known that if ^ = there 
is no tachyon mode, but there is a zero mode. 



19 



For < ^ < the potential still has the form of the left-hand panel of Fig. (volcano 
form). However, in this case the coefficient of the delta function is not large enough in order 
to support a tachyon mode. We can prove this result by solving numerically the eigenvalue 
equation (Eq. (34) in Ref. 0). We find that there is no tachyon mode in this region of ^. 

For ^ > the tachyon mode returns, as we obtain from Eq. (34) in Ref. . However we 
can not use Eq. (34) (of Ref. 151) if ^ > 1/5, as the index of the modified Bessel functions 
in Eq. (34)) becomes imaginary. In this case we have solved numerically the Schrondiger 
equation (31) by using a shooting method and we have confirmed that indeed the spectrum 
of the scalar field possesses at least one tachyon mode. In particular for ^ = + 1 we 
have exactly two tachyon modes, and in Fig. ^2 we have plotted the corresponding wave 
functions. Note that for ^ > the potential has the double-well form of the right-hand 
panel of Fig. ^2 In this region of ^ the spectrum of the scalar field is possible to contain 
one or more tachyon modes according to the depth of the double well. 

In this appendix we studied the spectrum of scalar field perturbations around the RS2- 
vacuum. Depending on the value of ^ we obtain that: (a) for ^ < we have unique tachyon 
mode, (b) for C > we have at least one tachyon mode and (c) for < ^ < there 
are no tachyon modes. We conclude that the RS2-vacuum is unstable against scalar field 
perturbations in cases (a) and (b). For ^ = or ^ = there are no tachyon modes, hence 
in these cases the RS-2 vacuum is stable. 

B Appendix: The tachyon character of the spectrum 
around the solutions of first and second class 

In this appendix we will show that the spectrum of scalar field perturbation around the 
solutions of first class (^ < 0) and second class (^ > ^c) exhibits a tachyon character. We 
can follow the analysis that is presented in the previous appendix. The only difference is 
that the value of the scalar field is nonzero in the bulk, hence the potential that appears in 
the Schrondiger equation (32) must be modified as 

V{w) = ^4^^ + a\w) (mw) + l2\<P\w)) (35) 

where w = Jq dz/a{z). 

In right-hand panel and left-hand panel of Fig. [T^l we have plotted the potential V{w) 
for the first class of solutions ^ < 0, and the second class of solutions ^ > respectively. As 
we see, the potential for both cases becomes negative and decreases without a lower bound. 
This implies a continuous spectrum of tachyon modes, and as a result both the first and 
second class of solutions are unstable. A possible way to resolve this instability is to put a 
second brane before the potential V{w) becomes negative. In this case, a fine tuning between 
the parameters of the model is necessary. 

C Appendix: Stable solutions for ^ = 

In this appendix we show that the solutions for C, = are stable against scalar field per- 
turbations. In the left-hand panel of Fig. Uni we have plotted the potential V{w) (see Eq. 
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Figure 12: The potential V{w)/k'^ as a function of kw^ for ^ = + 0.1 in the background 
of second class solutions (left-hand panel) and for ^ = —0.1 in the background of first class 
solutions (right-hand panel), for a/k = 6, A/k"^ = —6, ,A//c^ = 0.01. 

(35)) as a function of w. We observe that the potential is always positive and vanishes for 
a finite value of the coordinate w {kw = 1, 891). As we see in the left-hand panel of Fig. 
this finite value for w corresponds to infinite proper distance z. Thus, from the form of the 
potential V{w) in Fig. [121 we conclude that the scalar field spectrum consists of continuous 
modes with positive energies, which becomes discrete in the case of a second brane in the 
bulk. This result implies the stability of the solutions for ^ = ^c- 
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